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Models for structure formation attempt to predict the power spectrum of den- 
sity perturbations in the present universe from the initial power spectrum and 
the nature of dark matter. Observational constraints on the power spectrum 
CNj ■ at different scales in the present epoch can, therefore, be used to eliminate (or 
choose between) different theoretical models. Such a comparison is fairly easy 
! at large scales (at which linear theory is valid), and one can use observations 
like the MBR anisotropy, large scale steaming motions etc to constrain the 
Oh 1 models. But to discriminate between the models effectivley, it is necessay to 
constrain the power spectrum at small scales. The most reliable constraints on 
the power spectra at small scales come from the predicted abundance of bound 



systems which can be estimated reasonably accurately using Press-Schecter (or 
similar) methods 1 . In the past, this method has been used in conjunction with 
the quasar abundance 2-4 and cluster abundance 5-7 . We show here that the 
abundance of damped lyman alpha systems (DLAS, hereafter), provides a far 
stronger constraint on the models for structure formation. Models with a mix- 
ture of hot and cold dark matter 8-11 (which are consistent with large scale 
observations) are strongly ruled out by the DLAS constraints while models 
with cosmological constant 12 are marginally inconsistent. It is also possible 
to combine the constraints from the abundance of clusters, DLAS and QSO's 
to obtain model-independent bounds on the power spectrum at the nonlinear 
scales. These bounds are to be respected by any viable model for structure 
formation. 

The damped Lyman a systems have been studied extensively in the recent years. 
The spectroscopic surveys for DLAS have shown that 13-16 the absorbers - found between 
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redshifts 1.7 to 4 - have an average HI column density < N >~ 10 21 cm -2 . The most 
striking property of the DLAS comes from a determination of their abundance characterised 
by the mean number density of lines per unit redshift, dN /dz. Lanzetta et al (1991) find 16 
the number of damped Lya systems per unit redshift interval with N <; 2 x 10 20 cm -2 is 

0.16 ± 0.03 < dNdamp < 0,25 ± 0.04 at < z >= 2.5. (1) 

dz 

Here < z > refers to the average redshift of these systems and the two bounds are derived 
from two distinct subsamples of the data. Their data is also consistent with there being no 
evolution in the absorber properties. The number density of the damped systems exceeds 
that expected from present day galaxies by a factor of about 3 — 5 (see ref. 16). This 
has been interpreted to indicate that the cross section of the absorbers is larger than that 
associated with galaxies. 

More importantly, the observed value of dNdamp/ dz, used in conjunction with the 
mean column density allows one to estimate the density of neutral hydrogen contributed 
by the DLAS. The mean mass density of neutral gas contributed by the damped Lycn 
systems at their mean redhift < z > is pdam P (< z >) = [im v < N > (dN/cdt). Here 
\i = 1.4 is the mean molecular weight of the gas and m v is the proton mass. As the 
universe expands this average mass density would have decreased by a factor (1+ < z >) 3 . 
Comparing the resulting density with the present day critical density one gets the current 
density parameter of the HI making up the DLAS to be 

n damp = pm p <N> —\ <z> pc{1+<z>) [~ ) (2) 

(This result can be easily be modified to take into account a cosmological constant). To 
obtain the total mass associated with the DLAS we have to include two factors: Suppose 
/jv is the HI fraction of the gas in the absorbers and /& is the baryonic fraction in the 
universe. Then the density parameter contributed by the total mass associated with the 
DLAS, which has collapsed to form bound objects is 

q ^ ttdamp r 0-352 j f dN/dz \ f <N> \ (fNy 1 ( fb 

° LAS (f N fb) ~ (l + z) 1 / 2 V 0.2 J \W 21 cm- 2 J \0.5j VO-015/i" 2 J 

(3) 

for aO = l universe. At an average redshift of 2.5 one gets SIdlas = 0.094 if h = 0.5. 

It is not easy to estimate the mass associated with the individual DLAS directly. 
However, there are several indications to suggest that these systems have masses in the 
range of 1O 11_12 M . If we make this assumption that DLAS are galactic scale objects, then 
we see from (3) that their abundance will provide the most stringent constraint on models 
of structure formation : any such model should have enough power on galactic scales to 
be consistent with a collapsed fraction of about 10 percent, at the redshifts associated 
with DLAS. We elaborate on this point in greater detail below after considering briefly the 
evidence that the damped systems indeed represent galactic mass objects. 

For dN/dz = 0.2, a typical value indicated by (1) , one infers fldamp = 1.5/i _1 x 10~ 3 
for a flat universe. By comparison the mass density of the stars in present day galaxies 
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are thought 17 to account for fi„j s ~ 3/i -1 x 10 -3 . This similarity of 0^ amp and 0„i s first 
led to the suggestion that the damped Lycn systems may be the progenitors of present day 
luminous galaxies observed before the bulk of their stars formed, when the gas fraction was 
high. An alternate interpretation 18 is that these systems may represent a population of 
gas rich dwarf galaxies, which had smaller cross section but were much more abundant in 
the past. Direct information on the size and mass of the damped Lyct absorbers is limited 
to only a few cases at present: Study of the 21 cm absorption line 19 caused by the z = 2.04 
DLAS in PKS 0458-02, indicates that the absorber must extend more than 8/i _1 kpc across 
the line of sight. Spatial imaging 20 of the Lyct emission at z = 2.811 towards the QSO 
0528-250, shows that the DLAS at z = 2.811 could be a compact group of gas rich galaxies 
or a single large cloud of diameter ~ 100 kpc. Detection of CO emission lines 21 from the 
z = 2.14 DLAS in the same QSO indicates that the absorber may be a group of galaxies 
each galaxy having a gas mass ~ 2 x 10 n MQ and virial dimension ~ 40 kpc. Deep imaging 
of the field around the QSO 0000-263 has revealed a candidate galaxy associated with the 
z = 3.34 DLAS seen along the line of sight. 22 . In view of these observations, it appears 
likely that the damped Lya systems are gas rich massive galaxies rather than dwarfs. 
This is also consistent with the finding 23 that luminous galaxies, rather than dwarfs, are 
responsible for producing the Mgll absorption systems, of which the damped systems are 
a subset, at z = 0.3 — 0.9. We shall explore the consequences of this point of view below. 

To do this we have to compute the theoretically expected value of Q t heor y (M, z) con- 
tributed by bound systems with mass greater than M at any given redshift z. This can 
be done using the Press- Schecter formalism and we get the result that: 



^theory (> M,z) = erfc 



5c(l + z) 2 
V2a (M) 



(4) 



where erfc denotes the complementary error function, <jq(M) is the power spectrum eval- 
uated at z = by the linear theory and 5 C is the linearly extrapolated density contrast 
needed for collapse. In the simple spherical top hat model one finds that 5 C = 1.68. Given 
a theoretical model, we know <7 (M) and hence can estimate ^theory F° r a model to be 
viable we need Qdlas < ^theory, with the equality implying that every collapsed object 
is hosting a DLA system. This restriction, it turns out, can constrain the models very 
effectively. 

In figure 1 we have plotted O^ eory (lO 12 M0, z), the fraction of the mass which has 
collapsed into objects with 10 12 M© < M < 10 13 M©, for different theoretical models and 
compared them with the Qdlas(z) inferred from the data of Lanzetta et al. 16 , using (3) 
. The data points of Lanzetta et al. have been given by Turner and Ikuechi as a table 
of dN/dz at different redshifts. The data points with the error bars are for an O = 1 
universe. If O matter < 1, due to the presence of a non-zero cosmo logical constant, then the 
data needs to be shifted up. For 0„ ac = 0.8, the centers of the data need to be shifted to 
the location marked by the crosses. 

All the theoretical models in figure 1 are based on a primordial power spectrum of 
the form P(k) oc k. The solid curve is for the pure CDM model with the BBKS transfer 
function; the broken curve is for a model with CDM and cosmological constant (called 
V+CDM model herafter) with 0„ ac = 0.8, fl c dm = 0.2; the dotted curve is for the CDM 
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model with a bais factor of 2.5 (such a model, of course, cannot reproduce large scale 
observations including the COBE results unless the bias factor is scale dependent; we have 
included it only for the sake of comparison). The dot-dash curve is for a mixed H+CDM 
model with flhdm = 0.7, fi c d m = 0.3. All models (except, of course, the b=2.5 one) are 
normalised so as to reproduce the MBR quadrapole anisotropy measured by COBE. We 
have also taken h=0.5 for all cases except for the V+CDM model for which h=0.8. These 
parameters are chosen so that the models are consistent with the large scale observations. 

The comparison shows that the mixed dark matter models are strongly ruled out by 
the abundance of DLAS. The disparity is so high that the theoretical and observational 
uncertainties in the analysis is unlikely to be of any use in saving the model. The V+CDM 
model also falls short of producing the required abundance (note that in this case the data 
points are shifted up; it is the location of the crosses which one has to compare with the 
theoretical curves). However, the uncertainties in the data as well as some freedom which is 
available in theoretical modelling can be utilised to make the disparity lower. We feel that 
the data does rule this model also out though not as firmly and conclusively the H+CDM 
model. The pure CDM with COBE normalisation, of course, has no difficulty in explaining 
the abundances since it has a lot of small scale power; however this model is known to be 
inconsistent with other observations (like results from galaxy surveys). Finally, the (once 
popular) CDM model with b=2.5 also falls short of explaining the DLAS abundances. 

In the past, these models were tested as regards the observed quasar abundance. The 
constraints posed by the abundance of DLAS are far more stringent. To see this, let us 
calculate the Qquasan the density contributed by the mass associated with the quasars, 
along the folowing lines: The basic strategy is to relate the observed quasar abundance 
to the abundance of host galaxies of quasars. Quasars are thought to be powered by 
accretion on to black holes. Unless special models using super Eddington luminosities are 
involved, one can infer a characteristic black hole mass by assuming that the observed 
luminosity of the quasar corresponds to the Eddington luminosity. The mass of the host 
galaxy associated with a quasar of blue magnitude Mb can then be estimated to be: 

M G w ci x 1O 12 M x io°- 4 [- M s-26] (5) 

where 

i n/fbolx/ fb \—l,fhole\—l / r \ 

01 = L2( TfT ) W W } • (6) 

The three f-factors arise as follows: To convert the blue magnitude to the bolometric 
magnitude we have to use a correction factor, fool, which depends on the quasar spectrum 
but is in general of order of about 10. Secondly, only a fraction of matter in the universe 
may be in baryonic form. And, finally, only a fraction fhole of the baryons in a collapsing 
host galaxy be able to form the compact central object. We find from (5) that a bright 
quasar with Mb = —26 is typically hosted by a galaxy of mass of about 10 12 M Q . 
The density contributed by the quasars can now be written as 

rM B (M) 

n p „ snr (> M, z) = p- 1 / tM g (Mb)<S>(M b , z)dM B (7) 

J — oo 
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where Mb{M) is got by inverting (5) , $ is the quasar luminosity function and p is the 
mass density of the universe, at present. Further r = max(l,t(z)/tQ), takes into account 
the fact that only a fraction of order tq/t(z) of galaxies will display quasar activity if the 
quasar lifetime is smaller than than the age of the universe t(z) at redshift z. For the 
quasar luminosity function we shall use the form derived by Boyle et al. 25 for aO = 1, 
h = 0.5 universe. We shall also assume that their luminosity function, which is unevolving 
beyond z = 1.9, can be used upto z = 4 (see Ref. 4). Consider the abundance of collapsed 
halos required to explain the most luminous quasars, with Mb < —26. One can evaluate 
the integral in (7) numerically to find 

^auasar(> C^Mq, Z > 1.9) = 1.4 Cl X 10" 3 (8) 

For a particular model of structutre formation to successfully explain the abundance of 
quasars it must give ^theory (> M, z) > O gMasar (> M, z) . In figure 1 we have also plotted 
this line for 1.9 < z < 4.0. While this also rules out H+CDM model, it is clear that the 
constraint arising from this consideration is weaker than the one from DLAS. 

The abundance of objects like DLAS, quasars (and clusters) can be used to obtain 
bounds on the power spectrum ctq(M) at different scales. This is shown in figure 2 which 
is a plot of O versus a for different values of z. This plot is clearly independent of the 
power spectrum in the model and depends only on the validity of Press- Schecter analysis. 
The observed abundances of different objects at different redshifts can be used to mark out 
regions in this graph, thereby constraining the value of a at the mass scales corresponding 
to the objects. For example, the abundunce of Abell clusters at z ~ gives O ~ (0.001 — 
0.02); from the plot (marked by short horizontal lines in the z = curve) we see that this 
leads to the tight constraint of a ~ (0.5 — 0.7) at cluster scales [which turns out to be 
around 8/i -1 Mpc]. Similarly, the quasar abundance (marked by a straight horizontal line) 
gives a > (2 — 2.5) at galactic scales. The DLAS give a > (3 — 4) at the same scale (see 
the shaded region). Any model for structure formation should provide sufficient power at 
these galactic and cluster scales in order to be compatible with the observations. 

Similar constraints on the models for structure formation (from the abundance of 
DLAS) has been reached independently in ref. 26. 
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Figure Captions 

Figure 1: The density contributed by collapsed objects with mass in the range of 10 12 — 10 13 
in various theoretical models is compared with observations. All the theoretical models are 
based on a primordial power spectrum of the form P(k) oc k. The solid curve is for the pure 
CDM model with the BBKS transfer function; the broken curve is for a model with CDM 
and cosmological constant (called V+CDM model herafter) with Q, vac = 0.8, Q c dm = 0.2; 
the dotted curve is for the CDM model with a bais factor of 2.5. The dot-dash curve 
is for a mixed H+CDM model with flhdm = 0.7, O c ^ m = 0.3. All models (except, of 
course, the b=2.5 one) are normalised so as to reproduce the MBR quadrapole anisotropy 
measured by COBE. We have also taken h=0.5 for all cases except for the V+CDM model 
for which h=0.8. The data points with the error bars are based on the observed abundance 
of DLAS and are for an O = 1 universe. If O maWer < 1, due to the presence of a non-zero 
cosmological constant, then the data needs to be shifted up. For fl vac = 0.8, the centers 
of the data need to be shifted to the location marked by the crosses. The horizontal line 
is based on the abundance of quasars. 

Figure 2: The Q contributed by collapsed objects is plotted against the linearly extrap- 
olated density contrast <jq. The curves are parmetrised by the redshifts z = 0,1,2,3,4 
from top to bottom. The observed abundances of different objects at different redshifts 
are used to mark out regions in this graph, thereby constraining the value of a at the mass 
scales corresponding to the objects. The abundunce of Abell clusters at z ~ (marked by 
short horizontal lines in the z = curve) gives Q ~ (0.001 — 0.02) and leads to the tight 
constraint of a ~ (0.5 — 0.7) at cluster scales [which turns out to be around 8/i _1 Mpc]. 
Similarly, the quasar abundance (marked by a straight horizontal line) gives a > (2 — 2.5) 
at galactic scales. The DLAS give a > (3 — 4) at the same scale (see the shaded region) 
thereby providing the most stringent constraint. 
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